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Abstract — Unsteady interaction phenomena created hy the influence of the blade spacing have been reported in earlier 
experiments and CFD. However cyclic load changes in start-up and slow-down of the hydrodynamic torque converter 
operation have been beyond access to the current flow field calculation methods due to the extensive computer run time and 
memory requirements in the application of time dependent Navier- Stokes solvers to acceleration and deceleration. Therefore 
computations based on the use of ID mean-line flow simulation supported by optimised flow correction coefficients tuned in 
rig-test experiments enable to obtain solutions for engineering-type technical problems. 
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I. Introduction 


List of symbols 

Indices 

A 

cross section of the 

P pump 

E 

kinetic energy 

S stator 

c 

speed of the fluid 

T turbine 

CFD computational fluid dynamics 

1 st index number of blade row 

J 

moment of inertia 

2nd index inlet or outlet location 



1 inlet of blade row 

i 

swirl, angular momentum 

2 outlet of blade row 

M 

torque 


R 

radius vector 


V 

volume 


d 

outer diameter 


n 

revving speed 


t 

time 


il 

= pv = (M t n T ) / (M P n P ) = efficiency 


X 

conversion number 


P 

torque ratio = M out / M Pump 


V 

speed ratio n T / n P 


P 

density of the fluid 


(0 

angular velocity 
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II. Experience from the past 



Figure 1 View of the model W240 automotive hydraulic torque converter fabricated by 
Sachs showing the exit of the turbine blading on the outer diameter and the guide vane 

WHEEL DIRECTING THE FLOW FROM TURBINE TO PUMP NEAR THE HUB 

The behaviour of hydrodynamic torque converters has been reported by a number of investigators before. R. Herbertz [1] 
presents a non-linear system of partial differential equations to treat the changing speeds including acceleration and 
deceleration with the torque converter. E. Rohne [2] compares experiments involving unsteady torque under the presentation 
of steady-state conditions. He concludes that differences between time -dependent and steady-state results may occur at ratios 
v = n T / n P of revving speeds of turbine n T to pump n P which can originate from the influence of inertia caused by the rotating 
mass and an alternating circulation of fluid. In most investigations the moment of inertia of the rotating fluid is neglected 
because the unsteady phenomena are considered to be quasi-steady. Whenever the changes of the operating speed are slow, 
the quasi-steady operations fit closely to the description of the mathematically unsteady phenomena. The deceleration can be 
considered as rapid if 10 m/s 2 is exceeded. Values as high as 120 m/s 2 are attained during the sudden braking to stand-still 
from a large speed ratio to a small one with the fluid volume being 0.00285m 3 . 

Steady-state operation of hydraulic torque converters is explained and expressions are defined in VDI 2153 [3]. A set of 4 
dimensionless numbers (v, p, X, r|) characterises the behaviour of the physical operation from considerations of similarity in 
incompressible flow. The number p specifies the torque of the output shaft, driven by the turbine wheel, divided by the 
torque of the pump attached to the input shaft which drives the operating fluid by centrifugal force through the pump and into 
the turbine. The coupling of input and output shaft depends on the fluid flow momentum allowing for slip in the 
transmission. The torque conversion number p varies over the speed ratio v which is the ratio of turbine speed over pump 
speed. The power number X relates the torque of the pump with the 5th power of the diameter of the turbine rotor wheel the 
square of the angular speed and the fluid density: X = (M Pump ) / p co 2 Pump d 5 0u t P ut)- The efficiency r| is the product of the 

torque conversion number p and the speed ratio v representing the power ratio of turbine to pump r| = p v =(M T n T ) / (M P n P ). 
The brand name W240 originates from the German word "Wandler" for converter and the outer diameter of the turbine wheel 
of 240 mm. 

Early computer simulations in hydrodynamic torque converters are based on the ID mean-line flow theory as presented by T. 
Ishihara [4] and later by R. Leo [5]. The ID mean-line flow simulation relies on an ensemble of appropriately set coefficients 
representing total pressure loss and flow angle deviation which did not lead to acceptable results over the entire operating 
range unless an optimisation algorithm has been used for the data-point-wise adjustment of the ensemble of correction factors 
setting the difference between the computational result and the test rig result to zero over the entire operating range in 
repeated sweeps over corrections of deviation angle and total pressure loss. This optimisation is achieved with the algorithm 
of T. Rechenberg [6] in the computer code by H.P. Prufer [7]. 
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Figure 2 Four dimensionless numbers are v=speed ratio, p=power number, X=conversion number 
and T| EFFICIENCY CHARACTERISING THE BEHAVIOUR OF THE PHYSICAL OPERATION FROM 
CONSIDERATIONS OF SIMILARITY IN INCOMPRESSIBLE FLOW 


More refined extensions to unsteady computer formulations are presented by R. Herbertz [1]. Extensive efforts aimed at 
accessing hydrodynamic torque converter fluid dynamics with Navier Stokes solvers lead originally to the activity in SFB278 

[8] dealing with research on hydrodynamic power transmission cumulating finally in the investigation by M.Wollnik et al. 

[9] . 

Nevertheless, excessive computer run time requirements in unsteady 3D CFD led the way back to a ID approach for the 
problem of unsteady numerical simulation which gives an answer to the original question of SFB278 [8] in which E. Rohne 
[2] points out a hysteresis effect in the characteristics of the torque converter under steady-state considerations with alternate 
acceleration and deceleration. He concludes that differences in the characteristics may occur from the influence of inertia of 
the fluid sloshing about in the blading of the pump and the turbine wheels. Whenever the changes in operating speed are slow 
the quasi-steady computations fit closely to the description of the experiments in the test rig presented by R. Herbertz [1]. 
Investigations by P. De La Fuente [10] use the data basis of R. Herbertz [1] for unsteady operation as well as the data of H. J. 
Forster [11] for steady-state operation in a comparison. R. Herbertz [1] developes the system of equations for the unsteady 
torque conversion using a mean-line approach involving the rotational speed of the rotor rows and the mass flow balance 
from an unsteady momentum equation for each of the 3 cascade wheels and the energy conservation equation. In the case 
under investigation they are rotors of the pump, - of the turbine and of the freewheel, the latter with a speed range restricted 
to 0.8 < v < 1.0 in which the angle of incidence turns from positive to negative values [12]. 

Convective angular momentum 


— = Af. = — ffexc) pdV 
dt ‘ dt J. ’ 

Vi 

+ \(Rxd)pCdA 

Ai 

momentum in and out of the blade channel 


da) 

(lb) 


The sum of the angular momentum contributions of pump M P , turbine M T and guide vanes M Stator on all bladed wheels must 
go to zero: 


Mp + M t + M Statgr - 0 


( 2 ) 


Page | 169 


International Journal of Engineering Research & Science (IJOER) 


ISSN: [2395-6992] 


[Vol-2, Issue-4 April- 2016] 


The energy balance requires equilibrium between ingoing power in the pump, outgoing power from the turbine, dissipation 
by loss as well as the energy variation dE with changing speed over the time dt\ 


Pin + Pout + Ploss 


dE 

dt ( 3 ) 

Torque is influenced as well by the mass distribution incorporated in solid parts under rotation at varying speed of constant 
geometrical shape as by the more complicated case of the changing shape and speed of the volume of the work fluid 
transmitting the torque. The distribution of the fluid can be lumped into an appropriate shape representing the influence of 
inertia such as fluid in the shape of a rotating disc, ring or mass concentration positioned on representative radii. The moment 
of inertia J in the volume under rotation can be assumed to be of the shape of a disc with radius r and density p: 


/=J7 2 p</y 


M i = J i 


d co t 

dt 


( 4 ) 


Figures 3 and 7 show the situation with rapid acceleration, 8 and 9 with deceleration as the dimensionless torque rated by the 
torque M P0 and M T0 in the instant of the start at the initial moment t = 0. In Figures 7 the output shaft is unloaded. In this 
way the operating point goes from v = 0 and reaches v = 0.8 in 2.63 seconds. During the acceleration of v the motor remains 
with the pump at constant revving speed (1000 min' 1 ) at da> P / dt = 0. In Figure 11 the output shaft undergoes increased 
loading so that the turbine reaches stand-still in 0.23 seconds. Figures 3 and 1 1 show torque of the pump and the turbine as 
well as the speed ratio v versus time t. This allows determining the angular velocities oj, of pump and turbine by linear 
interpolation between the data points: 


dcoi 

rad 


f \ 

CO • .r\ CO • 

^ itL lt\ J 

> 1_ 

dt 

L 2 J 


A t [ 5 ] 




FIGURE 3 Measurement of rated torque M P /M P0 of 
the pump (- ♦ -) and of the turbine M T /M T0 (- 0 -) as 
well as the speed ratio v = n T / n P (-A-) at unsteady 
operation measured by R. Herbertz [1] at constant 
speed of the pump n P = 1000 min 1 under acceleration 
of the turbine during t = 2.63 s 



FIGURE 4 Comparison of rated turbine torque 
Mj/Mto in the unsteady experiment [1] with the 
steady-state computation. M T o is the initial torque 
from stand-still to acceleration (measurement - , 

computation -o-) 
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The ascend of the velocity ratio v is assumed to increase linearly between the data points taken during the measurement over 
time t (Figure 4). In this way it is possible to interpolate the angular speed at each instant between the respective data points. 
The acceleration changing with time may be described as the instantaneous angular velocity at the time the turbine requires 
from one data point to the next one. 

III. Quasi-unsteady treatment versus truly unsteady treatment 

In the following the comparison between the characteristics of the reference torque converter W240 with and without the 
influence of inertia is given as obtained from the transient experiment by R. Herbertz[l] including the moment of inertia as 
well as excluding it in the tests by H. J. Forster [11] using the steady-state computation. 

Figure 4 compares the computed steady-state case with the transient experiment during start-up in 3 seconds. The slower the 
variations are with time, the better is the agreement [2]. The case of acceleration shows satisfactory agreement between 
steady as well as unsteady treatment of the turbine. Considerable differences occur between prediction and experiment if the 
outlet shaft is decelerated to stand-still in 0.23s from v = 0.8 to v = 0. Figures 5a & 5b show a discrepancy between the 
unsteady experiment (- ♦ -) and the steady-state computer simulation (-o-) where the latter ignores the moment of inertia of 
rotating solid parts and fluid matter inside the bladed channel, thus increasing the torque of the turbine (figure 5b) and 
lowering the torque in the pump (figure 5a) with respect to the ID computer method for steady-state simulation. An extended 
theory following interaction of inertia of mass of solid and fluid matter follows the ID computational aproach with equations 
by R. Herbertz [1] and is assembled in a MATLAB programme [10]. 

IV. Unsteady computation using the mean-line theory 

The extension of the ID theory from the steady-state transmission of torque to the unsteady situations of acceleration, 
deceleration and change of load takes place by the interaction between inertia of rotating mass and fluid as well as the 
prescribed angular velocities of pump and turbine depending on the time as well as the retarding of pump and turbine wheel 
under the circumstances of the experimental results of the test rig. Using the optimised loss and deviation corrections with 
equations (6), (7) and (8) in table 1 below, minor deviations occur in the results of the characteristic curves (figure7). 


Table 1 

Momentum equations solved in the ID computer programme for pump (P), turbine (T), stator 

AND FREEWHEEL (S), ANY BLADED WHEEL (i), INERTIA OF SOLID BODY (J) AND FLUID (H) 
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Figure 5a 


Figure 5b 


FIGURE 5: Comparison of the rated torques M P /M P0 for the pump (above) and M T /M T0 for the turbine 
(below) obtained from unsteady measurements (- ♦ -) and steady-state computation (-o-) with a deceleration 

of the turbine in 0.23 s from v=0.8 to v=0 


The hydrodynamic torque is created by the fluid directed in the bladed channel in a way that angular momentum is 
transmitted. Furthermore the mass of the spinning rotors with the blading and pieces attached to input and output shafts 
contribute to the amount of torque acting on the shafts as soon as the revving speeds vary. A simple way to take into account 
the influence of inertia is to add torque from additional rotating angular momentum of solid parts M u and rotating fluid 
matter M m to consider the dynamic influence. 




(ii) 


The extension of the ID theory to allow for unsteady operation including acceleration, deceleration and change of load 
follows the approach of R. Herbertz [1] and is put down in a MATLAB computer programme solving the equations 6 , 7 and 
8 (table 1 below) under the boundary conditions of the time intervals as applied in the experiment making use of the 
characteristic operating lines of the optimised corrections for loss and deviation fitted to the results of the stationary 
experiments of the model W240 torque converter in the work by P. De la Fuente [10]. 



Figure 6 a 



Figure 5: Comparison of ID-computed torques of the turbine using unsteady (—A—) and the steady-state 
methods (-•-) for the turbine torque M T under an acceleration of 10 revs/s 2 (above) and of more than 100 
revs/s 2 (below) started from time t=0s. The torque of the pump changes less under the unsteady influence (- 

-A-, ~0~) 
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V. Situations under varying load 

Under a sudden increase of load on the turbine shaft the motor driving the pump reacts with a modest decrease of angular 
velocity. The deceleration of the pump remains at a rather constant -30.75m/s 2 (figure 6a). Figure 7 shows a result similar in 
trend to figure 6a with the torque in the unsteady as well as the steady condition of the accelerating turbine at 10 revs/s 2 but a 
discrepancy for the deceleration at 100 revs/s 2 starting at t=0 in figure 6b bas well as in figure 9. 

The example for the effect of inertia under acceleration shows a fairly similar trend of the rated turbine torque distribution 
M t /M to over the time, no matter how the computational prediction is done by either transient, steady-state, ID or 3D CFD 
computation as can be seen in figure 9 and the explanation in figure 14. Contrary to the acceleration in figure 9, the 
deceleration shows a wide spread in the distribution displayed in figure 10 for the distribution of rated torque over the time t. 
A rapid slow-down of the output shaft of the turbine causes the disruption of the unidirectional flow as the fluid speed 


overtakes the turbine wheel speed. 



Figure 7 Turbine torque M T /M T0 from transient 
measurements by R. Herbertz [1973] during acceleration 
from v = 0.01 to v= 0.8 in t- 2.63 seconds compared with 
the steady-state measurements by H. J. Forster [11]. 

- Transient- ( ) and steady-state measurement 

(-A-). 

- ID transient („ A „) and ID steady-state ( — A — ) 
computation and 

- 3D CFD computation transient ( ) and steady- 

state ( ) 

by P. De La Fuente [10] 

- Ratio of speeds n T /n P ( — A — ) 



FIGURE 8 Torque of pump (M P ) and turbine (Mr) 
measured over the time lapse of t-0.23 seconds during 
deceleration from v=0.08 to v=0.01 from unsteady 
measurements by R. Herbertz [1] and steady-state 
measurement by H. J. Forster [11] in comparison with 
the ID- and 3D-computer methods of P. De La Fuente 

[10] show significant discrepancies in torque M T /M T0 and 
Mp/Mpo during deceleration, in the measurements. 

Experiment: 

- Mp/Mpo transient ( ♦ ), M r /M ro transient (— o— ) , 

- Mp/Mpo steady-state (— n — ), M T /M T0 steady-state (— o-) , 

ID computation: 

- Mp/Mpo steady-state ( ), M T /M T0 steady-state ( ) 

- Mp/Mpo transient ( ), M T /M T0 transient (— o— ) , 

CFD computation: 

- Mp/Mpo steady-state ( ), M T /M T0 steady-state (— o-) 

- Mp/Mpo transient (— n — ), M T /M T0 transient (— o— ) , 
from P. De La Fuente [10] 
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f)l i 1 i — n- 

U 0 0.1 0.2 
t[s] 


FIGURE 9 Load applied to the output shaft by suddenly retarding the speed ratio from v=0.7 down to 
v=0.6 during 0.24 seconds obtained by an unsteady CFD result by P. De La Fuente [10] after a week of 
central processor time. The rated torque of the pump M P /M P 0 (-o-), - of the turbine M T /M T0 (-♦-) and the 
variation along the speed ratio v(-A-) are shown over the time t 


VI. Conclusion 

In the steady-state CFD analysis by M.Wollnik et al. [9] performance in computational effort required 1.4 Million elements, 
in the unsteady case. In the present 3D CFD simulation only 676 711 elements had been used for the simulation of a single 
blade channel in figure 8 and 9. Computer time of the order of one week central processor time had been required for the run 
indicated in figure 9 from v=7 down to v=6 during 0.24 seconds. Once the optimisation procedure has set the calibration 
coefficients correcting flow angle and pressure loss of the ID method in the curve fit to preceding rig tests, the ID 
performance prediction is applicable to the prediction of any torque converter of identical blade shapes and fluid mass at 
moderate computing times with better accuracy for the situation of acceleration than of deceleration (figures 7 and 8). 

At present the truely unsteady analysis on the basis of 3D Navier-Stokes solvers for the torque converter case does not seem 
economically justified for engineering applications because of excessive computer time requirement. 
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